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is  invoked  to  compute  macroscopic  stress-defonnation  responses  corresponding 
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choices  of  parameters  entering  a  macroscopic  plasticity  model,  with  different 
parameter  sets  corresponding  to  different  microstructures.  This  macroscopic 
model,  with  various  parameter  sets,  is  in  turn  used  to  simulate  impact  and 
perforation  of  a  thin  plate  of  the  aluminium  alloy  by  a  cylindrical  projectile. 
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1  Introduction 

Considered  here  is  the  dynamic  thermomechanical  response  of  a  monolithic  metallic 
plate  fabricated  from  a  precipitate-strengthened  aluminium  alloy.  The  alloy  of  study  is  A1 
2139,  consisting  primarily  of  elemental  A1  and  containing  smaller  amounts  of  Cu,  Mg, 
Ag,  Mn  and  trace  amounts  of  Si  and  Fe  (Cho  and  Bes,  2006).  The  heterogeneous 
polycrystalline  microstructure  exhibited  by  this  material  is  representative  of  the  general 
class  of  Al-Cu-Mg-Ag  alloys.  The  material  exhibits  a  face  centred  cubic  (FCC)  crystal 
structure  at  room  temperature  and  low  pressures,  and  consists  of  bulk  grains  and 
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secondary  phases  that  emerge  upon  processing  either  within  grains  or  preferentially  near 
grain  boundaries.  These  secondary  phases  may  include  Q  (chemical  composition  Al2Cu) 
that  nucleates  homogeneously  throughout  the  alloy,  contributes  to  hardening  without 
facilitating  grain  boundary  fracture  (Cho  and  Bes,  2006),  and  in  fine  distributions  leads  to 
creep  resistance  and  superior  strength  (Hutchinson  et  al.,  2001).  Also  in  such  alloys, 
without  the  addition  of  Ag  or  with  the  addition  of  Zr,  the  phase  denoted  O’  (also  of 
chemical  composition  Al2Cu)  emerges,  often  along  grain  boundaries  and  promoting 
intergranular  fracture.  Dispersoids  of  Mn  are  also  common.  Such  particles,  when  very 
small,  may  increase  ductility  by  promoting  mobile  dislocation  generation  (Vasudevan  et 
al.,  1986),  but  larger  Mn  particles  can  act  as  stress  risers  and  may  be  detrimental  (Brown 
and  Stobbs,  1971). 

Aluminium  alloys  of  various  compositions  have  been  the  focus  of  numerous 
numerical  studies  relevant  to  the  present  work.  Povrik  et  al.  (1995)  developed  a  crystal 
plasticity  model  incorporating  thermal  effects  to  describe  thin  film  aluminium 
interconnects.  Hanim  and  Klepaczko  (1999)  implemented  a  macroscopic  continuum 
theory  to  model  spall  fracture  in  an  aluminium  alloy.  Starink  et  al.  (1999)  described 
effects  of  precipitates,  ceramic  inclusions,  dislocations,  point  defects,  chemical 
compositions,  grain  sizes  and  subgrain  sizes  on  the  mechanical  behaviour  of 
Al-Li-Cu-Mg  alloys  using  a  scalar-based  phenomenological  model.  Gall  et  al.  (2000) 
studied  the  effects  of  particle  cluster  morphology  on  fracture  and  particle-matrix 
debonding  in  a  cyclically  loaded  cast  Al-Si  alloy.  Cheng  and  Itoh  (2004)  conducted 
simulations  of  ballistic  impact  of  thick  Al  7039  plates  using  macroscopic  plasticity  and 
failure  models.  Zhao  et  al.  (2004)  modelled  surface  roughening  in  aluminium  sheets 
under  both  uniaxial  and  biaxial  tension  using  a  computational  crystal  plasticity  approach. 
Dwivedi  et  al.  (2006)  modelled  shock  compression  of  polycrystalline  aluminium 
microstructures  featuring  a  variety  of  microstructural  features  such  as  hardened  or 
softened  grain  boundaries,  porosity  and  ceramic  inclusions.  Rezvanian  et  al.  (2006) 
studied  evolution  of  dislocation  cells  in  aluminium  undergoing  severe  plastic  deformation 
using  a  dislocation-based  crystal  plasticity  model.  Wu  et  al.  (2006)  modelled  grain 
boundary  damage  and  shear  localisation  across  multiple  grains  of  a  polycrystalline 
Al-Cu-Mg  alloy. 

A  limited  number  of  experiments  have  been  conducted  to  assess  ductility,  fracture 
strength  and  high-rate  behaviour  of  the  Al  2139  alloy  (Cho  and  Bes,  2006;  Ravichandran, 
2007).  Specifically,  the  alloy  exhibits  a  tensile  elongation  in  excess  of  0.15,  a  room 
temperature  fracture  toughness  on  the  order  of  200  MPa  m1/2,  and  an  ultimate  tensile 
strength  on  order  of  500  MPa.  At  high  compressive  strain  rates  (103  -  104/s),  the  material 
exhibits  relatively  little  rate  dependence  on  flow  stress,  with  a  maximum  flow  stress  on 
the  order  of  800  MPa  and  ductility  in  excess  of  0.5,  though  thermal  softening  is  evident  at 
compressive  strains  exceeding  0.35.  Observations  also  suggest  that  the  material  is 
relatively  resistant  to  failure  by  adiabatic  shear  localisation.  Preliminary  experiments 
(Cho  and  Bes,  2006;  Cheeseman,  2007)  have  indicated  promise  with  regards  to  resistance 
of  plates  of  the  alloy  to  dynamic  penetration. 

Crystal  plasticity  theory  is  used  in  the  present  approach  to  address  effects  of 
microstructure  -  specifically  rolling  textures  (Randle  and  Engler,  2000;  Wenk  and  Van 
Houtte,  2004),  slip  system  strengthening  rate  and  residual  elastic  energy  (Taylor  and 
Quinney,  1934;  Rosakis  et  al.,  2000)  on  the  dynamic  elastic-plastic  response  of  the  alloy. 
In  a  related  work,  Schoenfeld  and  Kad  (2002)  used  crystal  plasticity  theory  to  study 
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effects  of  crystallographic  texture  on  the  high-rate  shear  response  of  a  Ti-6A1-4V  alloy. 
Schoenfeld  and  Kad  (2002)  found  that  lattice  orientations  affect  the  overall  shear  strength 
of  the  metal  as  well  as  the  evolution  of  localisation  deformation  modes,  and  suggested 
that  the  most  shear-resistant  textures,  i.e.,  those  with  highest  macroscopic  shear  strength 
and  with  delayed  adiabatic  shear  localisation,  should  be  the  most  resistant  to  dynamic 
perforation.  Computational  complexity  presently  limits  3-D  crystal  plasticity  predictions 
to  domains  on  the  order  of  1  mm3  or  smaller  when  individual  grain  geometries  are 
fully  resolved  (Becker,  2004).  This  limitation  renders  simulations  of  impact  events 
with  specimens  of  realistic  spatial  dimensions  unfeasible.  As  such,  an  averaging 
procedure  (Taylor,  1938;  Kalidindi  et  al.,  1992)  is  used  here  to  deduce  the  effective 
thermomechanical  properties  of  an  aggregate  of  crystals,  and  then  these  results  are  used 
to  calibrate  macroscopic  model  parameters  entering  simulations  of  high-speed  impact.  Of 
course,  many  microstructural  details  are  lost  when  the  polycrystalline  geometry  is  not 
fully  resolved,  including  some  elements  of  anisotropy,  neighbouring  grain  interactions, 
shear  band  propagation  across  individual  grains  and  discrete  intra-  and  intergranular 
fractures.  But  such  effects  can  in  part  be  reflected  implicitly  via  choice  of  macroscopic 
thermodynamic  and  kinetic  material  parameters,  giving  what  is  deemed  here  a  reasonably 
realistic  depiction  of  the  polyciystalline  response  to  high-rate  shear  deformation.  The 
peripheral  shear  stress  along  the  lateral  impactor-target  interface  is  often  deemed  the 
dominant  resistive  stress  component  in  the  plugging  stages  of  perforation  (Schoenfeld 
and  Kad,  2002),  specifically  plug  formation,  plug  separation  and  plug  slipping  (Li  and 
Goldsmith,  1996).  It  is  noted  that  some  uncertainty  is  inevitably  inherent  in  the  present 
performance  predictions,  however,  as  a  result  of  the  limitations  imposed  by  the  relative 
simplicity  of  purely  macroscopic  depictions  of  inelasticity  and  failure. 

The  present  study  interfaces  aspects  of  materials  design  and  structural  analysis, 
making  efficient  use  of  limited  experimental  data  on  material  behaviour  and  dynamic 
performance.  The  two-scale  method  provides  quantitative  bounds  on  performance  that 
may  be  expected  via  conceivable  adjustments  of  the  alloy’s  microstructure.  The 
remainder  of  this  paper  is  organised  as  follows.  Section  2  discusses  the  single  crystal 
plasticity  formulation  used  for  microstructure  scale  modelling  of  precipitate-hardened  Al 
alloys.  Section  3  discusses  the  polycrystal  model  obtained  via  averaging  of  the  single 
crystal  response,  and  numerical  results  obtained  from  the  crystal  model.  Section  4 
presents  the  macroscopic  plasticity  fonnulation,  with  parameters  selected  following 
careful  consideration  of  the  aforementioned  crystal  plasticity  results.  Section  5  presents 
simulations  of  impact  and  perforation  of  the  Al  plate.  Conclusions  follow  in  Section  6. 

The  following  notational  scheme  is  used.  Scalars  and  individual  components  of 
vectors  and  tensors  are  written  in  italic  font,  while  vectors  and  tensors  are  written  in  bold 
font.  The  summation  convention  applies  on  repeated  indices,  with  3-D  Cartesian 
coordinates  implied.  The  •  symbol  denotes  the  scalar  product  of  vectors 

(a*b  =  aaba  =  axbx  +a2b2  +a363  j,  while  juxtaposition  of  second  rank  tensors  implies 
summation  over  one  set  of  adjacent  indices  j^(AB)“=  Aab  Bhc~^.  Summation  over  two  sets 

of  adjacent  indices  is  denoted  by  the  :  symbol  ( A  :  B  =  AahBab  j .  Superposed  •,  T  and  -1 

denote  the  material  time  derivative,  transpose  and  matrix  inversion  operations, 
respectively.  Additional  notation  is  explained  as  it  appears. 
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2  Single  crystal  plasticity 

The  crystal  plasticity  theory  implemented  here  follows  the  general  thermodynamic 
framework  developed  for  single  crystals  of  cubic  symmetry  by  Clayton  (2005a,  2005b), 
where  more  detailed  descriptions  and  physical  justifications  for  most  modelling 
assumptions  can  be  found.  The  framework  is  further  extended  here  to  encompass 
variability  in  properties  attributed  to  second -phase  inclusions  in  Al-Cu-Mg-Ag  alloys. 
The  second-phase  particles  of  interest  here,  on  the  order  of  tens  to  hundreds  of 
nanometres  in  diameter,  are  deemed  too  small  relative  to  the  dimensions  of  the  A1  grains, 
on  the  order  of  tens  to  hundreds  of  microns  in  size,  to  resolve  explicitly  and  distinctly 
from  the  bulk  grains.  Instead,  their  effects  are  assumed  to  enter  parameters  controlling 
kinetics  of  individual  slip  systems  (Rocks  et  al.,  1975)  and  mechanisms  of  energy  storage 
(Taylor  and  Quinney,  1934;  Rosakis  et  al.,  2000)  attributed  to  stress  concentrations  of 
blocked  dislocations  and  pile-ups  at  phase  boundaries,  for  example.  Because  precise 
values  of  such  parameters  have  yet  to  be  quantified  experimentally  for  the  present  class 
of  alloys,  parametric  studies  over  realistic  ranges  of  values  are  conducted,  ultimately 
enabling  description  of  effects  of  strain  hardening  rate  and  stored  energy  on  dynamic 
perfonnance  in  impact  events. 

Let  x  =  x(X,  t )  represent  a  smooth  time -dependent  motion,  with  x  and  X  denoting 
spatial  and  reference  coordinates,  respectively,  of  a  material  particle.  Deformation 
gradient  f  is  decomposed  multiplicatively  as: 

f  =  0x/0X  =  fefefp,  (1) 

where  f,  f0  and  f11  describe,  respectively,  kinematics  of  elasticity  and  rigid-body  rotation, 
stress-free  thermal  expansion  or  contraction,  and  the  cumulative  irreversible  contribution 
of  moving  crystal  defects.  In  general,  none  of  f®,  f8  or  f  is  integrable.  Elastic  and  thermal 
deformations  dictate  the  kinematics  of  the  slip  direction  vectors  stn)  and  slip  plane  normal 
covariant  vectors  nr“': 

s(a>  =  f 1 cf 1 es<“)  ,  m(c°  =  m^)fe_1fe_1 ,  (2) 

as  illustrated  in  Figure  1.  An  additive  split  of  spatial  velocity  gradient  I  follows  from  (1) 
as: 


1  =  Sx/<3x  =  ff-1  =  DT ^  +  f  ef  9f  +  f  ef  Bj  Pf  p-lf  8-lf  e-1 

=  le  Ti«  Tip 


(3) 


Thennal  deformation  of  a  cubic  lattice  is  isotropic: 

le  =  f  ef  s-1  =  aT0 1,  (4) 

where  9  is  the  temperature,  aT  is  the  temperature-dependent  thermal  expansion 
coefficient,  and  unit  tensor  is  denoted  by  1.  The  plastic  velocity  gradient  in  the 

intermediate  configuration  b  of  Figure  1  is  written: 


ip  =fpfp-i “)*{?) 
a= 1 


(5) 


120 


J.D.  Clayton 


with  7(a>  the  plastic  shearing  rate  on  slip  system  a,  n  the  number  of  potentially  active 
slip  systems,  and  <8>  the  tensor  product.  As  a  result  of  the  orthogonality  of  Sga)  and  , 
the  accumulated  plastic  deformation  fp  is  always  volume -preserving. 

Figure  1  Deformations,  configurations  and  slip  directors 


f  f° 


Localised  balance  laws  are  specified  in  the  current  configuration  as  follows: 

p  +  ptr{\)  =  0  ,  diva  +  pb  =  px,  a  =cT,  pe  =  o  :l-divq  + pr,  (6) 

with  p,  a,  b,  e,  q  and  r  denoting  the  current  mass  density,  Cauchy  stress  tensor,  body 
force  vector  per  unit  mass,  internal  energy  per  unit  mass,  heat  flux  vector  per  unit  current 
area  and  energy  source  per  unit  mass,  respectively.  Here  div  is  used  to  represent 
divergence  with  respect  to  spatial  coordinates,  and  tr  indicates  the  trace  of  a  second-rank 
tensor.  The  local  entropy  inequality  is: 

pf]>-div(ql  0)  +  pr  /  0,  (7) 

with  f]  the  rate  of  entropy  per  unit  mass.  Helmholtz  free  energy  per  unit  mass,  y/,  is: 

y/  =  e  -  0rj.  (8) 

Upon  substitution  of  (6)  and  (8)  into  (7),  the  entropy  inequality  becomes: 

o:\-{q‘Vx0)ld>p{y/  +  07i),  (9) 

with  V x0  the  spatial  temperature  gradient. 

A  Helmholtz  free  energy  potential  of  the  following  functional  form  is  assumed: 

¥  =  w(te,0^),  (10) 

where  the  intennediate  configuration  elastic  strain  ee  =(fe7fe  —  1) / 2.  The  symbol  £, 
denotes  a  dimensionless  scalar  internal  variable  representing  stored  micro-elastic  energy 
associated  with  crystal  defects  (e.g.,  sessile  dislocations)  that  may  impede  shearing  on 
each  slip  system.  Expanding  the  time  derivative  of  (10)  and  inserting  the  result  into  (9): 
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°  :  (le  + 1°  +  lp ) - (q-Vx0) /  0  >  p :  ee  +  {deyr) 8  +  {d^¥) i  +  0 77) ,  (11) 

where  the  subscript  following  the  5-operator  denotes  partial  differentiation  with  respect 
to  the  subscripted  variable.  Additional  algebraic  manipulations  give: 

^(se -p  +(aTtr(<s)~ p(dey/)- pfd 

(12) 

-P{d^+YJr(a)f(a)>^,0ie, 

a= 1 

where  p  is  the  mass  density  in  configuration  b,  and  the  elastic  second  Piola-Kirchhoff 
stress  se  =  y'efe_1afe_r  with  j e  =  detf' e  =  p  /  p.  The  resolved  Cauchy  stress  on  slip 
system  a  is  found  from  the  definition  T<a>  =  a  :  s<a>  ®  mla'.  The  following  relations  are 


deduced  from  (12): 

s  (13) 

t)  =  aTtr(a)l  p-dey/  =  %-dey/,  (14) 

(15) 

a= 1 

Rearranging  the  last  of  (6)  using  (13)  and  (14)  gives: 

n 

p0X~ p08e\f  =  pr-divq  +  'YJfa)f(a)  ~ p[dpp)S,.  (16) 

a= 1 

Specific  heat  capacity  c  is  defined  here  as  the  partial  derivative  of  the  internal  energy 
with  respect  to  temperature  at  constant  elastic  strain  and  constant  internal  variable,  i.e.: 

d  =  dee  =  -6dggy/  +  d^edgx,  (17) 


where  (8)  has  been  used.  Local  energy  balance  (16)  then  becomes: 

n 

pci)  =y~’r(a>x<a)  -  p[d pp  -  66  g£ip)j  +  p6d :ee  -  div  q  +  pr  .(18) 


temperature  a=\ 
change 


plastic 

dissipation 


energy  of 
lattice  defects 


thermoelastic  coupling 


heat 


heat 


conduction  supply 


A  free  energy  potential  per  unit  intermediate  configuration  volume,  deemed  applicable  to 
the  class  of  heterogeneous  alloys  of  study  here,  is  now  postulated  as: 


pW=X-te  :C:te  +  X-Kp[{\  +  0)^  +y. 


(19) 


where  S  =  tr(e?)  measures  the  elastic  volume  change,  C  is  the  temperature  dependent 
rank  four  elastic  modulus,  K\  and  p  are  temperature  dependent  elastic  stiffness 
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coefficients,  and  k  is  a  dimensionless  scalar.  The  scalar  Q,  also  dimensionless,  is  newly 
introduced  here  to  reflect  effects  of  composition,  e.g.,  solute  atoms  and/or 
second-phase  particles  in  aluminium  alloys,  on  stored  energy  in  the  bulk  crystal  lattice 
and  within  the  precipitates  themselves.  Function  y  = -pc9\n(9 1 90)  accounts  for 
thermal  energy,  with  90  a  reference  temperature  at  which  y  =  0.  From  (13)  and  (19),  the 
stress  satisfies: 

se=p3e,f=C:ec-^2l,  (20) 

where  the  second  term  on  the  right  of  (20)  accounts  for  materially  nonlinear  effects 
induced  by  elastic  volume  changes,  presumably  at  large  pressures.  The  hydrostatic 

pressure  p  and  deviatoric  stress  se'  in  the  intermediate  frame  then  become,  from  (20): 


fr(C:ee) 

3 


Kx92  ,se'  =se 


-  pi  =  C :  ee  - 


tr^l-Kx92  1.(21) 


Kinetics  of  plasticity  are  described  next.  A  power-law  viscoplastic  flow  rule  is  used  to 
address  deformations  induced  by  dislocation  glide  on  each  slip  system: 

(  ~(a)  V" 

r(a)=fo  ^j-  Sgn(f(a)).  (22) 

V& 

Here,  y0  and  m  are  material  parameters,  g(a)  is  the  slip  resistance  arising  from  dislocation 

barriers,  =  jeTi'a)  is  the  projected  shear  stress  pulled  back  to  the  intermediate 

configuration  b  of  Figure  1,  and  sgn(x)  =  x/|x|,  with  sgn(0)  =  1.  Thermal  softening 

attributed  to  increased  dislocation  mobility  at  high  temperatures  is  incorporated  by  the 
power  law  relation: 

g{a)  =  g(oa)  {o / 0O)P  ,  (23) 

with  the  flow  resistance  at  reference  temperature  0O  and  p  a  dimensionless  constant. 
The  following  relationship  between  the  ensemble-averaged  hardening  over  all  systems  at 
fixed  reference  temperature  and  the  internal  variable  if  is  used: 

~g(y'>  )  =  =  (24) 

a= 1 

with  an  initial  yield  stress,  b  the  magnitude  of  the  Burgers  vector  in  the  reference 

lattice,  and  pT  the  total  dislocation  line  length  per  unit  intermediate  configuration  volume 
associated  with  shearing  impedance.  Scalar  proportionality  factor  a  accounts  for 
dislocation  interactions.  Effects  of  intrinsic  lattice  resistance  to  dislocation  motion  (e.g., 
Peierls  barriers),  as  well  as  defects  present  initially  in  the  undeformed  lattice,  are 
incorporated  implicitly  in  the  value  of  g{‘“) . 

In  FCC  alloys  of  interest  here,  the  number  of  potentially  active  slip  systems  is  chosen 
as  n  =  12,  with  dislocation  glide  permitted  in  (110)  close-packed  directions  on  {111} 


Two-scale  modelling  of  effects  of  microstructure 


123 


planes.  Evolution  of  slip  resistance  at  reference  temperature  90  is  dictated  by  a  hardening- 
minus-dynamic-recovery  relation: 

g<,a>  =  A(\+0)^qap\j^\-Bg^Y\f^\  (25) 

p= i  p= i 

with  the  interaction  matrix  satisfying: 

qap=Sap+q{l-Sap),  (26) 


where  q  is  the  latent  hardening  ratio.  In  (25),  A  and  B  are  constants  and  <T>  is  a 
dimensionless  parameter  introduced  here  to  reflect  influences  of  alloy  composition  on  the 
slip  system  hardening  rates,  since  second-phase  particles  are  thought  to  drastically  affect 
strengthening  in  Al-Cu-Mg-Ag  alloys  (Hutchinson  et  al.,  2001;  Cho  and  Bes,  2006).  In 
metallic  crystals  containing  such  particles,  a  number  of  mechanisms  may  affect  the 
resistance  to  motion  of  individual  dislocation  lines  (Kocks  et  al.,  1975):  friction  stresses, 
anti-phase  boundaries,  misfit  stresses,  core  and  stacking  fault  energy  differences  and 
chemical  hardening.  Interfacial  mechanisms  at  particle-matrix  boundaries  may  lead  to  a 
flow  stress  dependence  proportional  to  pbcV2/w,  where  c  is  the  volume  fraction  of 
particles  of  width  w  (Kocks  et  al.,  1975).  In  (25),  phenomenological  scaling  factor  <T>  is 
applied  to  account  for  the  ensemble  effects  of  such  mechanisms  in  Al-Cu-Mg-Ag  alloys 
until  quantitative  experimental  data  becomes  available  to  suggest  a  more  descriptive 
formulation. 

In  the  numerical  implementation  of  the  model,  an  implicit  hyperelastic-viscoplastic 
algorithm  is  used  to  integrate  the  elastic-plastic  constitutive  response.  Slip  rates  for  a 
given  increment  spanning  times  t  and  t  +  At  are  found  implicitly  from  values  of  the 
resolved  shear  stress  and  hardening  variables  as: 


r{a)=h 


Aa) 


g%l 


(27) 


Notice  that  and  gj+l,  depend  upon  6  through  (23)  and  the  temperature  dependence 
of  elastic  moduli  in  (19).  In  an  adiabatic  analysis,  the  temperature  rate  for  a  given 
increment  spanning  t  and  t  +  At  is  found  from  (18),  in  the  absence  of  heat  conduction  or 
heat  sources,  as: 


pc  “ 
V  ^  a= 1 


(28) 


where 


P  = 


V 


YjT(a)f(a) -p(d4i//-Odg^i//)i  'YjT'-’y 


\a= 1 


(«).>(«) 


J\a= 1 


(29) 


such  that  1  /?  is  the  ratio  of  the  rate  of  stored  energy  to  plastic  dissipation. 

Additional  explanation  regarding  stored  energy  and  the  role  of  internal  state  variable 
£is  now  merited.  Manipulating  (24)  followed  by  time  differentiation  gives: 


124 


J.D.  Clayton 


~pd.il/ 


Kfi(l  +  ny  apn  ^ 

Furthennore  from  (19): 


""”^v  '  24  Pt  aPn7^\ 


g, oa>-  (30) 


~p9depi/ 


k(\+0)  6d0p  l  (a) 


ajLin 


k 

a= 1 


§6”'  -4a))- 


(31) 


Substituting  (30)  and  (31)  into  (29)  then  yields: 


(32) 


Finally,  the  dislocation  density  variable  can  be  determined  directly  from  the  first  of 
(30): 


P  =  1- 


/r(  1  +  /2)  ( /u  +  Ode/Li ) 


K*r-*r)z 


a=l 


a=l 


s 

a=l 


T(.a)-(a) 


Pr  - 


1 

apnb 


E(*r 


12 


(33) 


As  is  clear  from  (19)  and  (32),  positive  values  of  Q  amplify  stored  energy  associated 
with  internal  variable  £  and  equivalently,  residual  elastic  energy  engendered  by 
microscopic  stress  fields  manifesting  from  dislocation  density  pT.  Such  amplification 
would  manifest  physically  in  Al-Cu-Mg-Ag  alloys  from  the  presence  of  second-phase 
particles  of  Q,  O’  or  Mn  that  can  act  as  local  stress  risers  (Brown  and  Stobbs,  1971) 
or  barriers  to  dislocation  motion.  The  increase  in  stored  energy  would  be  attributed 
to  both  elastic  energy  of  the  deforming  second  phases,  misfit  stresses  (Kocks  et  al., 
1975),  as  well  as  elastic  self-energies  of  dislocations  that  may  build  up  at  phase 
boundaries. 

Parameters  entering  the  single  crystal  model,  many  specific  to  alloy  Al  2139,  are 
compiled  in  Table  1.  Single  crystalline  aluminium  exhibits  cubic  symmetry,  with  linear 
elastic  constants  Cn  =CU11,C12  =C1122,  and  C44  =C1212  obtained  from  Arsenlis  and 
Parks  (2002).  Derivatives  of  elastic  moduli  with  respect  to  temperature,  d&C,  are 
obtained  from  Flanim  and  Klepaczko  (1999).  Third-order  bulk  modulus  constant  Kh 
reference  mass  density  p0,  thermal  expansion  coefficient  aT,  melting  temperature  6U 
specific  heat  c  and  Burgers  vector  b  are  taken  from  Flatch  (1984).  The  value  of  the  latent 
hardening  ratio  q  is  obtained  from  Povrik  et  al.  (1995),  while  a  follows  from  Zhao  et  al. 
(2004).  Reference  shearing  rate  f0  and  strain  rate  sensitivity  m  are  obtained  from 
Rezvanian  et  al.  (2006),  with  the  relatively  large  value  of  the  latter  reflecting  little 
dependence  of  flow  stress  on  strain  rate.  Reference  temperature  90  is  chosen  for 
convenience  as  both  the  room  temperature  and  initial  temperature  during  adiabatic 
simulations.  In  fact,  all  simulations  discussed  later  in  Sections  3,  4  and  5  are  adiabatic, 
with  an  initial  temperature  of  300  K.  Adiabatic  heating  occurs  in  conjunction  with  plastic 
flow,  leading  to  well-known  thennal  softening  phenomena.  The  remaining  parameters  are 
obtained  via  calibration  to  experimental  data,  as  explained  in  Section  3 . 
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Table  1 

Parameters  for  crystal  model 

Parameter 

Value 

Parameter 

Value 

c„ 

108  GPa 

a 

2.3(1  Of5  +  1 .69(10)  8  A0/K 

Cn 

61.3  GPa 

c 

900  +  0.466  A0  J/kgK 

C44 

28.5  GPa 

sT 

155  MPa 

dgCU 

-6.9  MPa/K 

q 

1.0 

deCI2 

27.3  MPa/K 

a 

0.3 

dgC.44 

-17.1  MPa/K 

b 

0.2863  mu 

K, 

299  GPa 

fo 

0.001/s 

Po 

2700  kg/m3 

m 

50 

Ik 

10 

A 

480  MPa 

s0 

300  K 

B 

1.5 

933.5  K 

P 

-1.1 

The  procedure  used  here  provides  only  a  simple  link  between  precipitate  or  defect 
content,  stored  energy  of  the  lattice  and  yield  strength.  In  the  present  approach,  Q  and  0 
are  treated  as  adjustable  scalars  that  affect  energy  storage  and  slip  resistance  in  a  generic 
manner.  A  more  precise  approach  would  be  prescription  of  functional  fonns  for  Q  and  0 
for  the  material  that  depend  explicitly,  for  example,  upon  precipitate  or  dispersoid 
volume  fractions  and  sizes  (Kocks  et  al.,  1975),  internal  (residual)  lattice  straining 
associated  with  chemistry  and  interstitial  defects  (Mott  and  Nabarro,  1948),  and 
microscopic  debonding  and  fracture  mechanisms  associated  with  second  phases  and  their 
interfaces.  Such  a  more  precise  approach  will  be  enabled  when  sufficient  experimental 
data  characterising  the  microstructure  and  relevant  microscopic  behaviour  of  Al  2139 
becomes  available. 


3  Polycrystal  plasticity  and  homogenised  response 

The  crystal  model  is  calibrated  to  available  high-rate  compression  data  and  then  used  to 
predict  the  high-rate  shear  response  important  in  plugging  stages  of  perforation  (Li  and 
Goldsmith,  1996;  Schoenfeld  and  Kad,  2002).  A  computational  algorithm  is  used  to 
generate  average  stress-strain  and  temperature  histories  for  aggregates  of  many  grains, 
first  in  order  to  calibrate  unknown  constitutive  model  parameters  specific  to  alloy  Al 
2139,  and  then  to  predict  such  histories  for  feasible  variations  in  parameters  associated 
with  different  processing  routes.  In  these  deformation-controlled  simulations,  constant 
velocities  are  prescribed  on  boundaries  of  a  cube  of  material  of  unit  dimensions,  with  this 
cube  containing  one  or  more  grains  of  the  material.  Taylor’s  assumption  (Taylor,  1938; 
Kalidindi  et  al.,  1992)  is  invoked  for  polycrystalline  simulations,  meaning  that  f  of  (1)  is 
identical  for  each  grain  in  the  aggregate  throughout  the  time  history  of  deformation.  High 
strain  rates  are  considered  here,  such  that  adiabatic  conditions  apply  and  local  thermal 
interactions  such  as  heat  conduction  between  crystals  are  not  addressed.  The  average 
Cauchy  stress  o  is  calculated  for  the  aggregate  for  each  time  increment  in  the  analysis 
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via  volume  averaging,  assuming  each  crystal  occupies  an  equal  part  of  the  total  aggregate 
volume  V  in  the  reference  configuration: 

o  =  V~ljadV.  (34) 

Figure  2  shows  the  model  predictions  for  an  aggregate  of  300  initially  randomly  oriented 
crystals,  deformed  under  uniaxial  stress  conditions  at  constant  compressive  strain  rates  of 
£■  =  3500/8  and  £  =  8 100/s.  In  coincident  reference  and  spatial  coordinates,  the  total 

deformation  gradient  is  prescribed  by  f  \  =  1  -  st  where  e  is  a  constant  and  the 

logarithmic  strain  is  e  =  -  In  /[ .  Also  shown  in  Figure  2  is  the  smoothed  experimental 
stress-strain  response  obtained  from  a  compressive  split  Flopkinson  bar  experiment 
(Ravichandran,  2007)  at  a  rate  of  s  =  8100/s.  The  compressive  stress  in  Figure  2  is  the 

only  nonzero  component  of  the  average  stress  of  (34),  and  is  labelled  as  a  =  a1 1 .  In  the 
simulations  of  Figure  2,  grains  in  the  material  are  specified  as  randomly  oriented  since 
after  final  heat  treatments,  quenching  and  aging,  A1  2139  plates  are  not  substantially  cold 
worked  or  rolled  to  yield  a  preferential  texture  (Cho  and  Bes,  2006),  though  experimental 
texture  analysis  would  provide  a  stronger  justification  for  the  modelling  assumption  of 
initially  random  orientations.  Three  hundred  was  found  to  be  a  sufficient  number  of 
grains  to  provide  an  isotropic  stress-strain  prediction  that  did  not  vary  substantially  with 
further  increases  in  the  number  of  crystals. 

Figure  2  Average  compressive  stress  versus  logarithmic  compressive  strain:  experiment  and 
polycrystal  model 


In  high-rate  experiments  (Ravichandran,  2007)  variations  in  strain  rate  from  ~107s  to 
~104/s  produced  little  changes  in  flow  stress,  implying  that  the  deviatoric  strength  of  A1 
2139  is  relatively  rate  independent  over  this  range  of  strain  rates.  The  rate  independence 
is  manifested  in  the  crystal  plasticity  model  via  the  large  value  of  m  entering  (22)  and 
Table  1,  and  is  demonstrated  by  the  close  agreement  in  stress  predictions  at  rates  of 
£•  =  3500/s  and  s  =  8100/s  in  Figure  2.  Plasticity  model  parameters  g^,  A,  B  andp  are 

selected  in  order  to  obtain  the  consistency  between  model  and  experiment  demonstrated 
in  Figure  2,  with  <P  of  (25)  taken  as  zero  by  default.  Notice  that  the  material  exhibits  a 
peak  compressive  strength  on  the  order  of  850  MPa  at  a  logarithmic  compressive  strain  of 
approximately  0.35.  Subsequently,  stress  relaxation  occurs,  emerging  in  the  model  via 
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both  dynamic  recovery  in  (25)  and  thermal  softening  in  (23).  The  latter  phenomenon 
occurs  as  plastic  work  is  converted  to  temperature  rise  in  adiabatic  analyses  via  (28). 


Figure  3  Projectile-plate  geometry  in  simple  plugging  scenario,  with  rolling  coordinate  system 


{  impact  velocity 


Table  2 

Ideal  textures  for  rolled  FCC  crystals 

Texture 

{ hkl}(uvw ) 

Euler  angles  \i/,0,<j> 

Cubic 

{001}  (100) 

0,  0,0 

Copper 

{ii2}(nT) 

180,35,  135 

S 

{123}  (634) 

211,37,  117 

Brass 

{no}  (T12) 

35,  90,  135 

Goss 

{110}(001) 

0,  90,  135 

{100} 

{001}  (ran) 

Random,  0,  0 

The  crystal  model  is  next  used  to  predict  the  response  to  high-rate  shearing  deformation, 
specifically  resistance  to  shearing  in  the  direction  tangential  to  the  lateral  impactor-plate 
interface  occurring  during  the  plug  formation  and  slippage  portions  of  perforation.  As 
shown  in  Figure  3,  this  shear  stress  component  cr13  is  labelled  simply  as  r.  Also  shown  in 
Figure  3  is  the  coordinate  system  typically  used  to  describe  rolled  textures  (Randle  and 
Engler,  2000),  with  RD  denoting  the  rolling  direction  X1,  TD  denoting  the  transverse 
direction X2  and  ND  denoting  the  normal  direction  X3.  In  addition  to  the  ensemble  of  300 
randomly  oriented  grains  described  already,  five  single  crystalline  orientations 
corresponding  to  textures  obtained  during  conventional  rolling  operations  on  plates  or 
sheets  of  FCC  metals  were  investigated:  cubic,  copper,  S,  brass  and  Goss  textures.  The 
corresponding  crystallographic  orientations,  specifically  the  Euler  angles  ( S)  in 
Roe’s  notation,  are  listed  in  Table  2,  following  Wenk  and  Van  Floutte  (2004).  In  Table  2, 
{hkl}  denotes  the  crystallographic  plane  normal  aligned  in  the  normal  direction  ND,  and 
(mvw)  denotes  the  crystallographic  direction  aligned  parallel  to  the  rolling  direction  RD. 
The  sixth  texture  listed  in  Table  2,  labelled  as  {100},  consists  of  300  grains  with  one 
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cube  axis  aligned  parallel  to  the  ND,  but  with  the  other  two  cube  axes  randomly  rotated 
about  the  ND.  This  provides  for  transversely  isotropic  symmetry  in  the  {100} 
microstructure,  in  contrast  to  the  other  more  conventional  rolling  textures  that  generally 
exhibit  orthorhombic  symmetry  (Wenk  and  Van  Houtte,  2004). 

Figure  4  shows  adiabatic  shear  stress-strain  behaviour  for  the  randomly  oriented  and 
textured  A1  2139  crystals.  Shearing  is  conducted  in  the  one  to  three  plane,  i.e., 
/  3  =  1  +  ft.  The  macroscopic  accumulated  shear  strain  is  then  found  simply  as  y  =  ft, 

while  the  average  shear  stress  f  =  V~x  j" cr13dV.  The  prescribed  strain  rate  of  f  =  104/s  is 

a  magnitude  typically  observed  during  dynamic  perforation  (Schoenfeld  and  Kad,  2002; 
Clayton,  2005a).  From  Figure  4,  crystals  with  cubic  orientation  are  predicted  to  offer 
substantially  greater  resistance  to  shear  deformation  than  the  other  textures  at  strains 
y  <  0.6,  while  conversely,  crystals  with  copper  orientation  are  predicted  to  offer 
substantially  less  resistance  to  shear  deformation  at  strains  y  <  0.3.  The  aggregate  with 
{100}  texture  exhibits  slightly  less  shear  strength  than  the  cubic  texture,  but  is  still 
substantially  stiffer  than  the  remaining  cases  for  y  <  0.6.  It  should  be  noted  that  the 
predictions  for  various  textures  are  not  obtained  via  calibration  to  experimental  data  on 
single  crystal  orientations  since  such  data  remains  unavailable  for  A1  2139;  hence,  some 
uncertainty  exists  in  the  predictions  of  Figure  4.  Instead,  as  explained  above,  the  model 
parameters  are  calibrated  to  polycrystal  data  assuming  random  orientations  (Figure  2), 
and  differences  in  response  among  different  textures  then  arise  naturally  in  the 
computations  as  a  result  of  elastic  and  plastic  anisotropy,  the  latter  induced  by  initial  and 
evolving  orientations  of  the  slip  directions. 

Figure  4  Average  shear  stress  versus  shear  deformation  at  f  =  104/s:  polycrystal  model  with 
initial  textures 


In  what  follows  next,  thermomechanical  behaviours  of  eight  different  microstructures, 
corresponding  to  eight  sets  of  crystal  parameters  -  specifically  variations  in  texture, 
hardening  rate  and  residual  energy  storage  -  are  investigated.  Microstructures  and 
parameters  are  listed  in  Table  3.  Microstructures  1  and  5-8  feature  random  textures  with 
300  grains,  while  microstructures  2  and  3  feature  cubic  and  copper  textures,  respectively, 
corresponding  to  the  stiffest  and  most  compliant  orientations  with  respect  to  high-rate 
shearing  as  indicated  in  Figure  4.  Microstructure  4,  an  aggregate  of  300  grains  randomly 
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oriented  about  the  ND  axis  ({100}  orientation),  is  also  studied.  Microstructures  5  and  6 
include  nonzero  values  of  0,  with  microstructure  5  corresponding  to  a  50%  increase  in 
the  rate  of  hardening  associated  with  the  first  term  on  the  right  of  (25),  and  microstructure 
6  corresponding  to  a  50%  reduction  in  this  rate  of  hardening.  Differences  in  hardening 
rates  may  be  attributed  to  various  concentrations  of  alloying  ingredients  in  Al-Cu-Mg 
alloys  leading  to  different  solutions  and  precipitate  distributions,  for  example,  as 
discussed  by  Starink  et  al.  (1999),  Hutchinson  et  al.  (2001)  and  Cho  and  Bes  (2006).  Note 
that  the  material  is  not  considered  free  of  precipitates  or  associated  lattice  defects  when  0 
=  0.  The  value  0  =  0  is  used  as  only  as  a  relative  reference,  with  enhancement  of  strain 
hardening  due  to  such  defects  represented  by  0  >  0,  and  reduction  in  hardening 
represented  by  0  <  0.  Finally,  microstructures  7  and  8  exhibit,  through  nonzero  values  of 
Q,  increased  energy  storage  capacity.  As  will  be  shown  later,  positive  values  of  Q  result 
in  smaller  values  of  the  heat  dissipation  fraction  /?  of  (29)  and  (32),  leading  to  less 
adiabatic  heating  than  would  occur  for  otherwise  equivalent  microstructures  with  Q  =  0. 
As  mentioned  already,  stress  concentrations  at  phase  boundaries  as  well  as  residual 
stresses  supported  directly  by  second-phase  particles  in  Al-Cu-Mg-Ag  alloys  presumably 
may  contribute  to  variations  in  Q. 


Table  3  Microstructures  investigated  in  polycrystalline  and  macroscopic  simulations 


Microstructure 

Texture 

0 

Q 

1 

Random 

0.0 

0.0 

2 

Cubic 

0.0 

0.0 

3 

Copper 

0.0 

0.0 

4 

{100} 

0.0 

0.0 

5 

Random 

1.5 

0.0 

6 

Random 

0.5 

0.0 

7 

Random 

0.0 

1.0 

8 

Random 

0.0 

2.0 

Average  effective  stress-shear  strain  behaviours  for  polycrystalline  aggregates 
corresponding  to  the  eight  cases  of  Table  3  are  shown  in  Figure  5.  Again  the  prescribed 
deformation  is  pure  shear,  f\  =l  +  jt,  at  a  constant  rate  of  ^  =  104/s.  The  effective 
Cauchy  stress  follows  the  usual  definition: 


with  individual  stress  components  on  the  right  side  of  (35)  following  from  volume 
integral  (34).  Note  that  (35)  is  consistent  with  the  notation  used  in  Figure  2,  since  the 
effective  stress  and  the  normal  stress  are  equal  under  the  uniaxial  conditions  considered 
in  Figure  2.  Note  also  that  the  shear  stress  definition  used  in  Figure  4  and  the  effective 
stress  of  definition  (35)  are  related  by  <r  =  v/3r.  From  Figure  5,  the  following 
observations  are  made.  Microstructures  2  and  4  exhibit  the  largest  shear  strengths  (i.e., 
flow  stresses)  at  strains  y  <  0.4  due  to  respective  cubic  and  {100}  orientations,  and 
microstructure  5  exhibits  the  largest  shear  strength  at  strains  y  >  0.4  due  to  the  prescribed 
increased  hardening  rate.  Microstructure  3  exhibits  the  lowest  shear  strength  among  all 
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cases  at  strains  y  <  0.17,  because  of  its  copper  texture,  while  microstructure  6  exhibits  the 
lowest  shear  strength  at  larger  strains  y  >  0.17  because  of  its  prescribed  reduced 
hardening  rate.  Microstructures  7  and  8  demonstrate  slightly  higher  flow  stresses  than  the 
base  microstructure  1,  because  of  the  reductions  in  temperature  rise  and  correspondingly 
reduced  thermal  softening  in  the  former.  Recall  that  in  micro  structures  7  and  8,  increased 
energy  storage  leads  to  less  dissipation  and  less  adiabatic  temperature  rise.  For 
microstructure  7,  the  increase  in  flow  stress  relative  to  microstructure  1  begins  to  emerge 
at  y  ~  0.2,  while  such  increases  emerge  slightly  sooner  in  the  deformation  history,  at 
strains  around  y  ~  0. 17,  for  microstructure  8. 

Figure  5  Average  effective  stress  versus  shear  deformation  at  y  =  104/s:  polycrystal  model  with 
various  microstructural  parameters 


Figure  6  Average  ratio  of  net  dissipated  energy  to  plastic  work  versus  shear  deformation  at 
f  =  104/s:  polycrystal  model  with  various  microstructural  parameters 


Thermodynamics  of  energy  storage  are  investigated  in  more  detail  in  Figure  6. 
Considered  in  Figure  6  are  microstructures  1,  7  and  8  of  Table  3,  corresponding  to 
randomly  oriented  polycrystals  with  values  of  Q  of  0,  1  and  2,  respectively.  The  results 
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correspond  to  the  same  high-rate  simulations  described  in  Figure  5.  The  average  heat 
dissipation  fraction  shown  in  Figure  6  is  computed  from: 

P  =  V~X\  fdV,  (36) 

where  the  local  value  of  f  for  each  crystal  of  the  aggregate  follows  from  (32). 

For  microstructure  1  with  Q  =  0,  a  range  of  0.8  <  <1.0  over  the  strain  history  is 

obtained,  a  range  typically  observed  and/or  assumed  for  many  FCC  metals  undergoing 
significant  plastic  deformation  (Taylor  and  Quinney,  1934;  Aravas  et  al.,  1990).  Such  a 
result  is  used  to  motivate  the  choice  of  the  value  of  k  listed  in  Table  1.  For  micro  structure 
7  with  Q  =  1,  the  lower  bound  of  f  drops  relative  to  that  of  microstructure  1  such  that 
0.55  <  /?<1.0  over  the  deformation  history  in  the  former.  For  microstructure  8  with 
Q=  2,  the  lower  bound  drops  even  further  and  0.02  <  /?  <  1.0.  In  each  case,  the  predicted 
value  of  f  decreases  initially  with  strain  as  dislocations  accumulate  in  conjunction  with 
hardening  mechanisms,  leading  to  an  increased  rate  of  energy  storage  associated  with 
self-equilibrated  internal  stress  fields  supported  by  the  defects.  Then  at  larger  strains 
y  >  0.3,  an  increase  in  f  with  strain  is  observed  as  the  rate  of  dislocation  accumulation 
slows  due  to  dynamic  recovery  prescribed  in  (25).  Similar  trends  and  ranges  of  f  falling 
between  bounds  predicted  here  for  Q  =  1  and  Q  —  2,  have  been  observed  in  rate- 
insensitive  alloy  Al  2024-T351  (Rosakis  et  al.,  2000). 


4  Macroscopic  model  formulation  and  calibration 

A  macroscopic  thermoviscoplastic  model  is  implemented  for  use  in  impact  simulations  of 
a  steel  cylinder  striking  a  monolithic  flat  plate  of  the  Al  2139  alloy.  The  macroscopic 
plasticity  model  parameters  for  the  Al  plate  are  calibrated  repeatedly  to  stress-strain 
results  from  the  polycrystal  model  shown  in  Figure  5,  with  a  different  set  of  parameters 
obtained  for  each  microstructural  feature  that  is  varied  among  simulations  listed  in 
Table  3:  hardening  rate,  energy  storage  rate,  or  texture.  The  macroscopic  damage  model 
for  the  Al  plate  is  calibrated  to  experimental  observations  of  tensile  and  compressive 
failure  for  use  in  most  of  the  impact  simulations.  However,  effects  of  variable  ductility  of 
the  Al  plate  on  its  dynamic  performance  are  also  studied  via  systematic  adjustment  of  the 
damage  model’s  parameters.  The  macroscopic  framework  also  relies  upon  a 
Mie-Gruniesen  equation  of  state  (Johnson  et  al.  1997)  providing  the  relationship  between 
volumetric  strain,  pressure  and  temperature  and  accounting  for  elastic  nonlinearity  at  high 
pressures.  In  a  certain  sense,  the  deviatoric  plasticity  component  of  the  present  work  can 
be  considered  multiscale  since  macroscopic  parameters  emerge  from  finer  scale,  crystal 
plasticity  modelling,  while  the  volumetric  response  and  damage  model  are  purely 
macroscopic  in  origin. 

The  viscoplastic  response  follows  from  an  associated  flow  rule  whereby  the  radius  of 
the  isotropic  yield  surface,  in  tenns  of  scalar  deviatoric  stress  measure  (35),  is  given  by: 

a  =  [Y0+YlENp][\-(9iex)M], 


(37) 
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with  Yq  the  initial  yield  stress,  ep  the  cumulative  equivalent  plastic  strain,  and  other 
nonzero  material  constants  Yu  N  and  M  reflecting  linear  and  exponential  strain  hardening 
and  thermal  softening,  respectively.  Fundamentally  different  forms  of  flow  rules  have 
been  suggested  in  (22)-(23)  and  (37).  Either  type  of  rule  is  sufficient  to  reproduce  the 
limited  experimental  data.  The  microscopic  flow  rule  of  Section  2  was  developed  for 
anisotropic  single  crystals,  while  the  macroscopic  rule  of  Section  4  is  appropriate  only  for 
isotropic  polycrystals.  The  simple  macroscopic  rule  (37)  is  chosen  here  for  computational 
efficiency,  without  stringent  regard  for  the  physics  of  plastic  deformation.  More  complex 
models  with  a  stronger  basis  in  physics  are  available  (Wright,  2002),  but  the  limited 
amount  of  experimental  data  that  exist  for  A1  2139  is  insufficient  to  enable  unique 
selection  of  all  parameters  entering  more  detailed  models. 

The  three-dimensional  implementation  features  an  additive  decomposition  of  the 
spatial  deformation  rate  into  elastic  and  deviatoric  plastic  parts,  and  explicit  numerical 
integration  using  a  radial  return  algorithm  with  the  Jaumann  rate  of  Cauchy  stress  taken 
as  the  objective  rate.  Additional  details  regarding  the  numerical  implementation,  mostly 
standard  in  the  computational  plasticity  literature,  can  be  found  in  Johnson  et  al.  (1997). 
In  the  implementation  of  the  macroscopic  model,  energy  equation  (18)  is  simplified  as: 

pcd  =  /?o  :lp  +  g,  (38) 

where  lp  is  the  spatial  plastic  velocity  gradient  (or  the  plastic  strain  rate  since  o  =  o7),  /? 
is  a  constant,  and  the  function  f  accounts  for  temperature  changes  manifesting  from  the 
equation  of  state. 

In  simulations  of  impact  events  involving  large  deformations,  failure  models  are 
typically  needed  to  capture  the  loss  of  strength  commensurate  with  shear  localisation, 
fracture,  or  spall.  In  the  model  implemented  here,  damage  is  assumed  to  accumulate  in 
the  material  with  cumulative  inelastic  deformation,  with  the  maximum  tolerable  amount 
of  damage  accumulation  increasing  with  the  stress  triaxiality;  i.e.,  damage  accumulation 
is  accelerated  by  tensile  pressure  and  retarded  by  compressive  pressure.  Upon  failure,  the 
mass  of  the  material  is  conserved,  but  the  solid  is  presumed  to  no  longer  support 
deviatoric  stress  or  tensile  hydrostatic  stress.  The  cumulative  inelastic  strain  sP  at  which 
failure  occurs  is: 

4  =Dof(0)exp(DlP/a),  (39) 

where  D0  is  the  failure  strain  under  conditions  of  pure  shear,  and  D\  accounts  for 
increases  (decreases)  in  ductility  with  increasing  (decreasing)  pressure  p  =  -tr(a)/3. 
Finally,  /  accounts  for  possible  effects  of  temperature  on  ductility.  In  the  numerical 
implementation,  (39)  is  used  in  an  incremental  fonn  since  the  local  pressure  and 
deviatoric  stress  state  may  vary  with  the  time  history  of  deformation.  Note  that  more 
elaborate  and  physically  descriptive  plasticity  and  damage  models  could  lend  more 
realism  to  the  macroscopic  simulations  when  sufficient  experimental  data  become 
available  to  populate  the  extensive  parameter  sets  of  more  complex  models.  Because 
sufficient  data  are  not  available  to  calibrate  /  for  Al  2139,  here  /  =  1  is  used  in  all 
simulations,  regardless  of  microstructure. 

In  the  impact  simulations  that  follow,  values  of  parameters  entering  (37)-(39)  are 
varied  in  order  to  reflect  microstructures  modelled  explicitly  with  crystal  plasticity  theory 
in  Section  3.  As  listed  in  Table  4,  values  of  parameters  entering  (37)  are  varied  in 
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calibrations  of  the  macroscopic  plasticity  model  to  the  high-rate  results  of  Figure  5, 
assuming  a  stress  state  of  pure  shear.  Recall  that  the  material  was  found  to  be  fairly  rate 
insensitive  over  the  range  of  high  strain  rates  of  present  interest;  hence  the  absence  of 
strain  rate  effects  on  hardening  in  (37)  is  thought  justified.  In  simulations  1-8,  constants 
corresponding  to  weighted  averages  of  f  introduced  in  (36)  over  the  deformation  history 
are  used  in  (38),  with  emphasis  given  to  values  at  strains  prior  to  failure  y  <  0.3.  Values 
of  Y0,  Y],  N  and  M  are  then  chosen  such  that  the  results  of  (37)  and  (38)  best  fit  the 
polycrystal  plasticity  predictions  of  Figure  5. 


Table  4  Simulations,  microstructures  and  corresponding  macroscopic  model  parameters 


Sim.  # 

Microstructure 

Yo(MPa) 

YfMPa) 

N 

M 

P 

A 

A 

1 

1 

610 

2,250 

0.23 

1.10 

0.945 

0.3 

2.0 

2 

2 

300 

3,000 

0.10 

1.30 

0.954 

0.3 

2.0 

3 

3 

300 

3,500 

0.40 

1.20 

0.947 

0.3 

2.0 

4 

4 

250 

3,000 

0.15 

1.40 

0.952 

0.3 

2.0 

5 

5 

100 

4,000 

0.30 

1.40 

0.838 

0.3 

2.0 

6 

6 

1,630 

125 

0.20 

0.80 

0.999 

0.3 

2.0 

7 

7 

610 

2,250 

0.26 

1.23 

0.790 

0.3 

2.0 

8 

8 

610 

2,600 

0.30 

1.20 

0.552 

0.3 

2.0 

9 

1 

610 

2,250 

0.23 

1.10 

0.945 

0.45 

2.0 

10 

1 

610 

2,250 

0.23 

1.10 

0.945 

0.15 

2.0 

For  microstructures  featuring  randomly  oriented  grains,  use  of  an  isotropic  material 
model  is  deemed  acceptable,  at  least  for  monotonic  (as  opposed  to  cyclic)  defonnation 
scenarios  where  directional  back-stresses  need  not  be  delineated.  However,  for 
microstructures  2-4,  the  shearing  response  of  plastically  anisotropic  crystals  is  used  to 
calibrate  an  isotropic  macroscopic  model.  Such  an  assumption  may  over-  or 
under-estimate  the  deviatoric  strength  of  the  material  in  other  directions.  On  the  other 
hand,  if  the  shearing  response  dominates  resistance  to  perforation  and  plugging  as  has 
been  suggested  elsewhere  (Li  and  Goldsmith,  1996;  Schoenfeld  and  Kad,  2002),  the  use 
of  an  isotropic  macroscopic  model  to  compute  the  resistance  of  textured  sheets  to 
penetration  may  at  least  provide  reasonable  bounds  on  performance  that  suggest 
directions  for  material  improvements,  if  not  precise  results.  More  realistic  results  for 
textured  sheets  could  presumably  be  attained  by  calibration  of  the  crystal  predictions  to 
an  anisotropic  macroscopic  theory,  for  example  one  featuring  an  orthotropic  yield  surface 
(Hill,  1949).  Finally,  for  cylindrical  projectiles  of  the  sort  considered  here,  shear  stresses 
resisting  penetration  may  be  non-uniform  along  the  circumference  of  the  projectile-target 
interface  for  rolling  textures  lacking  transverse  isotropy.  Hence,  the  results  reported  for 
the  sheet  with  transversely  isotropic  { 100}  texture  (simulation  4)  are  thought  to  provide  a 
more  realistic  prediction  of  perfonnance  for  a  rolled  sheet  with  cubic  texture  than  those 
of  simulation  2  that  assume  all  shear  stresses  occur  tangential  to  {001}  planes  in  the 
lattice. 

Remaining  material  parameters  are  selected  as  follows.  In  simulations  1-8,  values  of 
D0  and  D\  are  chosen  based  on  experimental  observations  (Cho  and  Bes,  2006; 
Ravichandran,  2007)  that  the  A1  2139  exhibits  a  tensile  elongation  of  -0.15,  a  shear 
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ductility  of  ~0.3,  and  a  compressive  ductility  in  excess  of  0.5.  In  simulations  9  and  10, 
respectively,  the  ductility  of  the  material  is  increased  or  decreased  by  a  factor  of  0.5,  in 
order  to  quantify  effects  of  ductility  on  simulated  performance.  Library  parameters 
typical  for  aluminium  alloys  (Kohn,  1969)  are  used  for  the  macroscopic  equation  of  state 
of  A1  2139,  details  of  which  are  available  elsewhere  (Johnson  et  al.,  1997).  Possible 
affects  of  variations  of  microstructure  on  the  hydrostatic  response  are  not  addressed. 

The  above  discussion  pertains  to  the  modelling  of  the  target  Al  plate.  The  steel 
projectile’s  behaviour  is  modelled  identically  in  all  simulations,  using  standard 
appropriate  versions  of  (37)-(39)  accounting  for  high  strain  rate  effects  on  plasticity  and 
fracture,  an  equation  of  state  for  the  volumetric  response,  and  available  macroscopic 
library  parameters  for  4340  steel  (Kohn,  1969;  Johnson  and  Cook,  1983,  1985).  Also  for 
the  steel,  the  default  /?  =  1  is  used  in  (39). 


§  Impact  simulations 

The  simulations  considered  here  involve  a  projectile  comprised  of  4340  steel,  impacting  a 
square  Al  2139  plate  of  19.05  mm  (0.75”)  thickness.  The  projectile  is  modelled  as  a  right 
circular  cylinder  (Gama  et  al.,  2001),  of  diameter  20  mm,  length  22  mm  and  mass  53.7  g. 
The  lateral  dimensions  of  the  plate  are  an  order  of  magnitude  greater  than  its  thickness, 
i.e.,  190.5  mm  x  190.5  mm.  The  large  lateral  dimensions  ensure  that  the  perforation 
process  concludes  before  stress  wave  reflections  from  the  boundaries  affect  the  behaviour 
of  the  projectile.  The  initial  velocity  of  the  projectile  in  all  simulations  is  specified  as 
1,000  m/s,  perpendicular  to  the  surface  of  the  target  (i.e.,  null  obliquity).  The 
performance  metric  chosen  here  is  the  residual  velocity  vR  of  the  projectile  after 
perforation:  the  smaller  the  residual  velocity,  the  better  the  performance  of  the  plate.  The 
present  approach,  in  which  the  plate  is  severely  overmatched  by  the  projectile,  renders 
quantitative  performance  appraisal  more  straightforward,  via  consideration  of  Vr,  than 
would  be  the  case  if  perforation  did  not  occur  in  all  instances.  This  approach  also  reduces 
possible  dependence  of  the  solution  on  mesh  density  that  may  become  excessive  as  the 
limit  velocity  approaches  zero.  In  order  to  obtain  a  precise,  mesh-independent  prediction 
of  the  limit  velocity,  a  constitutive  and  numerical  formulation  including  an  intrinsic 
length  scale  of  dimension  larger  than  the  grid  spacing  should  be  used,  for  example  one 
resolving  micro  structural  features  or  incorporating  spatial  gradients  of  field  variables 
such  as  plastic  deformation  or  damage  (Voyiadjis  and  Abu  Al-Rub,  2006). 

Several  mesh  densities  were  considered,  with  only  half  of  the  full  geometry 
represented  as  permitted  by  symmetry  reductions.  The  performance  metric  was  found 
to  be  nearly  insensitive  to  the  number  of  nodes  and  elements  in  the  mesh  at  the  chosen 
mesh  density,  consisting  of  a  total  of  138,240  tetrahedral  elements.  A  dynamic 
Lagrangian  analysis  package  was  used,  with  default  values  for  parameters  controlling 
time  integration,  artificial  viscosity  and  element-to-particle  conversion  (Johnson  et  al. 
1997).  Regarding  the  latter,  in  order  to  obviate  problems  encountered  with  large  element 
distortions,  finite  elements  undergoing  very  large  deformations,  i.e.,  effective  strains 
greater  than  0.5,  are  converted,  via  a  mass-conservative  transformation,  to  a  meshless 
environment  with  a  generalised  particle  algorithm  (Johnson  et  al.,  2000). 

Results  from  simulation  1  of  Table  4  are  visualised  in  Figure  7,  at  four  time  instances 
of  the  impact  event.  The  steel  cylinder  clearly  undergoes  large  compressive  deformation 
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along  its  axial  direction,  and  lateral  expansion  in  the  vicinity  of  the  impact  face.  Plugging 
failure  of  the  aluminium  plate  is  also  observed,  with  the  velocity  of  the  plug  clearly 
exceeding  that  of  the  projectile  at  t  =  300  /us  in  Figure  7(d).  In  addition  to  the  primary 
plug,  a  debris  cloud  of  smaller  fragments  of  the  A1  alloy  is  also  visible  behind  the  plate. 
The  results  depicted  in  Figure  7  for  simulation  1  are  qualitatively  similar  to  those 
observed  in  the  other  cases  listed  in  Table  4. 

Figure  7  Simulation  of  perforation  of  microstructure  1  at  (a)  t  =  0  (b)  t  =  20  /us  (c)  t  =  100  /us 
(d)  t  =  300  /is 

VVA 

(a)  (b)  (c)  (d) 

Figure  8  Average  final  residual  velocities  of  projectile  for  ten  simulations  corresponding  to  ten 
different  hypothetical  microstructures  of  the  A1  2139  plate 
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Residual  velocities  for  all  nine  simulations  are  shown  in  Figure  8.  These  are  obtained 
from  the  solution  data  via  the  spatial  integral: 

Vr  =V~1fv3ciV’  (40) 

with  V  the  instantaneous  (i.e.,  current  configuration)  volume  of  the  projectile  and  v3  the 
component  of  the  velocity  vector  in  the  normal  direction  (ND)  to  the  plate,  as  introduced 
in  Figure  3.  Residual  velocities  remained  relatively  constant  at  simulation  times  t  >  1300 
/us,  though  the  values  shown  in  Figure  8  correspond  to  the  instant  t  =  300  /us.  Recall  that 
each  simulation  corresponds  to  a  different  set  of  material  parameters  (Table  4)  associated 
with  variations  in  micro  structure  (Table  3).  The  baseline  result  for  comparison  with  the 
other  cases  is  simulation  1,  featuring  isotropic  properties  corresponding  to  a  random  set 
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of  grain  orientations,  with  strength  and  ductility  calibrated  most  closely  to  available 
experimental  data  (Cho  and  Bes,  2006;  Ravichandran,  2007).  For  simulation  1, 
vR  =  525  m/s. 

In  general,  the  greater  the  strength  or  ductility  of  the  plate,  the  more  energy  absorbed 
by  the  plate  during  impact  and  perforation,  the  greater  the  reduction  in  kinetic  energy  of 
the  cylindrical  projectile,  and  lower  the  cylinder’s  residual  velocity.  Simulations  2  and  4 
exhibit  very  low  residual  velocities,  vR  =  258  m/s  and  vR  =  300  m/s,  respectively,  as  a 
result  of  increased  shear  strengths  attributed  to  texturing.  These  correspond  to 
microstructures  2  and  4  of  Figure  5. 

Because  of  its  high  hardening  rate,  simulation  5  also  exhibits  a  lower  residual  velocity 
than  case  \,vR  =  495  m/s,  though  the  reduction  in  residual  velocity  is  not  as  drastic  as  for 
simulations  2  and  4  because  the  largest  increases  in  strength  in  simulation  5  manifest  only 
at  large  strains,  visible  in  Figure  5  for  microstructure  5  with  y  >  0.4.  According  to  (39) 
and  the  failure  model  parameters  used  here,  such  large  strains  are  not  sustainable  by  the 
material  unless  the  stress  state  is  highly  compressive.  In  the  present  simulations,  the 
material  tends  to  lose  strength  and  energy  absorbing  capacity  prior  to  realisation  of  the 
increased  strength  that  would  be  afforded  at  very  high  pressures. 

Simulations  3  and  6  suffer  degradation  in  performance  relative  to  base  simulation  1 
because  of  their  reduced  shear  strengths,  translating  to  a  reduction  in  kinetic  energy 
mitigation  of  the  plate.  Recall  from  Figure  5  that  in  simulation  3  the  strength  is  reduced 
due  to  lattice  orientation  effects,  while  in  simulation  6  the  strength  is  reduced  due 
to  a  decreased  hardening  rate  associated  with  composition,  e.g.,  solution  effects  or 
second-phase  particle  content  and  morphology.  As  such,  the  residual  velocities  are 
relatively  high:  in  simulation  3,  vR  =  637  m/s  and  in  simulation  6,  vR  =  599  m/s. 

Simulations  7  and  8  offer  the  possibility  of  improved  dynamic  performance  of  the  A1 
2139  alloy  via  an  increase  in  the  material’s  ability  to  store  residual  elastic  energy.  In 
these  cases,  increases  in  residual  elastic  energy  result  in  a  decrease  in  the  adiabatic 
temperature  rise,  leading  to  increased  strength  since  thermal  softening  and  shear 
localisation  processes  (Wright,  2002)  are  impeded.  However,  in  the  present  investigation, 
simulation  7  yields  a  comparable  residual  velocity  to  simulation  1,  since  benefits  in 
strength  attributed  to  increased  elastic  energy  storage  are  small  and  do  not  emerge  in 
simulation  7  until  relatively  large  shear  strains  are  attained.  Simulation  8,  on  the  other 
hand,  with  greater  energy  storage  (Figure  6)  and  strength  increases  at  lower  strains 
(Figure  5)  provides  modest  improvement  in  performance  relative  to  case  1,  with  vR  =  486 
m/s. 

Finally,  effects  of  ductility  are  quantified  by  comparing  results  from  simulations  1,  9 
and  10.  A  50%  increase  in  ductility,  assuming  no  change  in  the  influence  of  triaxiality 
upon  failure  strain,  provides  a  substantial  reduction  in  residual  velocity  and  a  noticeable 
improvement  in  performance,  with  vR  =  409  m/s  in  simulation  9.  Conversely,  a  50% 
decrease  in  ductility  results  in  an  expected  increase  in  residual  velocity  and  inferior 
performance,  with  vR  =  648  m/s  for  simulation  10. 


6  Conclusions 


A  two-scale  approach  has  been  implemented  to  estimate  effects  of  fine  scale  material 
parameters  on  the  macroscopic  constitutive  behaviour  of  an  Al-Cu-Mg-Ag  alloy  and  its 
resistance  to  impact  and  perforation.  At  the  fine  scale  of  resolution,  a  crystal  plasticity 
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model  has  been  developed,  accounting  for  finite  deformations,  grain  orientation  effects 
and  consistent  thermodynamics  associated  with  temperature  changes,  lattice  defects  and 
residual  elastic  energy.  Flexibility  has  been  inserted  into  the  model  enabling  description 
of  the  effects  of  composition  on  slip  system-level  hardening  and  energy  storage,  for 
example,  attributed  to  interactions  between  dislocations  and  second-phase  precipitates. 
An  averaging  scheme  has  been  used  to  obtain  macroscopic  model  parameters 
corresponding  to  hypothetical  microstructures  of  Al-Cu-Mg-Ag  alloys  conceivably 
attainable  via  processing  operations.  Such  operations  include  texturing  by  rolling  or 
strengthening  via  addition  of  alloying  ingredients.  Macroscopic  simulations  of  impact 
and  perforation  of  thin  plates  of  the  alloy  were  conducted,  and  effects  upon  performance 
of  material  properties  corresponding  to  different  crystalline  microstructures  were 
quantified  in  terms  of  residual  velocities  of  the  projectile. 

The  results  suggest  that  the  most  substantial  improvements  in  dynamic  perfonnance 
of  the  A1  alloy  may  be  achieved  by  rolling  to  instil  a  cubic  texture  or  by  taking  processing 
steps  to  substantially  increase  the  alloy’s  ductility,  e.g.,  maximum  elongation  increases 
on  the  order  of  50%.  Additional  benefits  may  be  achieved  from  composition  changes  that 
provide  increases  in  rates  of  slip  system  hardening,  e.g.,  an  increase  in  strengthening  rate 
on  the  order  of  50%,  or  increases  in  the  rate  energy  storage,  e.g.,  a  reduction  in  the 
average  fraction  of  dissipated  plastic  work  from  0.95  to  0.55.  In  the  present 
theoretical-computational  study,  each  microstructural  variation  has  been  examined  in 
isolation,  though  property  variations  are  likely  highly  coupled  in  physically  realistic 
structure -processing  paths.  For  example,  adjustment  of  strength  of  the  alloy  without 
alteration  of  its  ductility  may  prove  difficult  using  traditional  processing  routes. 
Flowever,  the  present  research  provides  some  quantification  of  improvements  to  impact 
and  perforation  resistance  that  could  be  realised  via  realistic  adjustments  to  the 
microstructure  of  A1  2139,  and  some  guidance  as  to  which  microstructure  properties  and 
features  may  be  of  highest  priority  when  considering  the  dynamic  performance  of 
metallic  plates  from  a  micromechanics  perspective.  Clearly  the  most  substantial  benefits 
would  emerge  from  processing  efforts  that  seek  to  simultaneously  increase  strength, 
ductility  and  energy  storage  mechanisms. 

It  is  noted  that  a  number  of  limitations  lend  uncertainty  to  perfonnance  predictions  of 
the  crystal  plasticity  model  and  subsequent  macroscopic  impact  simulations.  The  crystal 
plasticity  predictions  should  be  viewed  as  most  accurate  only  for  loading  conditions  for 
which  the  model  has  been  calibrated,  namely  polycrystals  defonned  in  compression, 
initially  at  room  temperature  at  strain  rates  up  to  approximately  104/s.  The  robustness  of 
the  crystal  plasticity  model,  and  the  quality  of  the  calibration  of  parameters,  could  be 
improved  via  consideration  of  additional  experimental  data  obtained  for  crystals 
subjected  to  various  stress  states,  and  in  various  non-random  initial  orientations. 
Regarding  the  ballistic  simulations,  experimental  data  used  for  calibration  of  deviatoric 
plasticity  parameters  and  the  macroscopic  failure  model  are  limited  to  macroscopic  data 
collected  over  a  range  of  strain  rates  up  to  104/s,  even  though  higher  strain  rates  may  be 
achieved  in  the  target  material  during  the  impact  event.  The  macroscopic  models  used 
here  do  not  attempt  to  explicitly  resolve  the  discrete  nature  of  failure  mechanisms  such  as 
fracture  and  shear  localisation  that  could  be  addressed  more  realistically  in  a  multiscale 
context  (Clayton,  2006)  if  accompanied  by  sufficient  experimental  characterisation.  For 
three  of  the  cases  investigated  here  featuring  textured  materials,  some  effects  of 
anisotropy  are  neglected  in  the  macroscopic  model,  specifically  the  stress  components 
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other  than  that  component  providing  the  shear  resistance  to  plugging  on  planes 
orthogonal  to  the  rolling  direction.  Hence,  the  benefit  of  cubic  texturing  (simulation  2) 
reported  in  Figure  8  may  be  overestimated,  with  the  transversely  isotropic  {100}  texture 
(simulation  4)  offering  perhaps  a  more  reasonable  prediction  of  the  maximum 
performance  improvement  attainable  via  rolling  of  the  A1  plate.  Nonetheless,  the  present 
results  do  provide  previously  unavailable  guidance  as  to  which  micromechanical 
phenomena  should  be  given  high  priority  when  developing  alloys  for  resistance  to  high 
velocity  impact  and  perforation. 
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